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ABSTRACT 



Modifications were made to Corgnati's (1971) steady state, symmet- 
rical hurricane model, which was based on a model proposed by Riehl 
(1963), and in which bulk aerodynamic transfer formulas were used to 
predict the air-sea interaction below a specified vortex. The modifi- 
cations resulted in a more flexible model in which hurricane intensity 
was controlled by sensible and latent heat fluxes, which were in turn 
dependent on the ocean surface temperature. A boundary layer model 
based on similarity theory was then used, to specify air-sea interaction, 
by replacing the bulk transfer formulas with Cardone’s (1969) extension 
of Blackadar’s (1965) two-layer neutral, baroclinic boundary layer 
model. To achieve realistic results the time scale was modified to 
include centripetal acceleration, which resulted in a decreasing height 
of frictional influence with decreasing radius. Varying Cardone’s for- 
mulation of roughness length and the ratio of eddy exchange coefficients 
for heat and moisture resulted in stress and heat flux values in agree- 
ment with values determined from hurricane observations. Including the 
complete boundary layer model resulted in greater sensitivity to ocean 
temperature variations than when bulk transfer formulas were used to 
determine the air-sea interaction. 
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I. INTRODUCTION 



It is well recognized that sea surface temperature affects hurricane 
intensity; the hurricane winds in turn, exert a stress on the ocean sur- 
face and induce the extraction of heat and moisture. A "steady state" 
hurricane model requires that a balance exists in this energy exchange 
between air and ocean. 

Riehl (1963) described a two-layer model for an axisymmetric hurri- 
cane in which a balance existed between the heat source and momentum 
sink at the air-sea interface. The "steady state" upper level vortex 
was then specified based on the air-sea interaction within the boundary 
layer. 

Corgnati (1971) used Riehl 1 s model to construct a numerical model 
which used predictive equations for heat and moisture within the 
boundary layer, and bulk transfer formulas for fluxes at the air-sea 
interface. It was felt that considerable improvement in Corgnati* s 
model could be attained through the use of recent boundary layer models 
based on Monin-Obuhkov "similarity theory"; specifically through the 
use of Cardone’s (1969) extension of Blackadar's (1963) two-layer neutral 
baroclinic boundary layer model. 

The purpose of this study was to modify Corgnati* s numerical model 
so as to better define the air-sea interaction. 

Several numerical simulations of hurricane development have been 
made with complicated multi-level models. These models include some 
interaction with the ocean. The scope of this study was to model only 
the boundary layer beneath an upper level vortex, defined by Riehl *s 
model, which responds to the air-sea interaction. 
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II. DISCUSSION OF THE MODEL 



A. REVIEW OF RIEHL' S MODEL 



Riehl' s two-layer model consisted of an inflow layer and an outflow 
layer (Figure 1). Riehl and Malkus (1961) have shown the conservation 
of absolute angular momentum is a reasonable first approximation in the 
outflow layer. This constraint leads to the expression for computing 
the radius at which the tangential component of the outflow goes from 
cyclonic to anticyclonic. 



R o = (2 V V f) 

1 



1/2 



( 1 ) 



where R. is the eye wall radius and V is the tangential wind at R.. 

i 

Riehl and Malkus (1961) also demonstrated the approximate conserva- 
tion of potential vorticity in the inflow layer. Riehl (1963) applied 
this to the equation of motion in polar coordinates (see Figure 2) 



dV 



V V 
e r 



dt + R 



+ fV = - 



1 dP 



p rS 0 + p bz 
and integrated over the depth of the inflow layer to get 

RT = constant , 



1 © 



( 2 ) 



(3) 



where is the tangential surface stress. Using the expression 

T e = C D P s V e /cos > (4) 

s 

and assuming the drag coefficient (C^) and surface density (p g ) to be 
constant, also that the inflow angle is small enough so that cosine (X 
can be closely approximated by unit yields the expression 

R^^ - constant . (5) 
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Figure 1. Illustrating hurricane mass circulation V q is tangential 
wind component, 0^ is equivalent potential temperature, R is radius. 




Figure 2. Illustrating polar coordinate system used, V 0 is tangential 
wind component, V is radial wind component, R is radius, 9 is angular 
measurement. 
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Riehl (1963) and Shea (1972) have shown this to be a reasonably good 
approximation to the V q -R profiles. 

Equation (2) can also be solved for the radial wind profile, 



V 

r 



AP 




( 6 ) 



where G is the acceleration of gravity, and f is the coriolis parameter. 

Integration of the hydrostatic equation, for typical values of 
equivalent potential temperature (0^) in a hurricane, produces the 
relationship 

P' = - 2.56 9^ , (7) 



where P’ is the pressure departure from 1005 mb and 0^ is the departure 

from 350K. Riehl (1963) also determined a relation between the maximum 

tangential wind (V ) and the equivalent potential temperature (0 ) at 

°i e i 

the eye wall by equating V to the height gradient of the isobaric 

i 

surfaces and applying a cyclostrophic relationship to obtain 

V e (m/sec) = 14.1(6 - 6 ) 1/2 , (8) 

i e i e o 

where subscript "o" denoted the value at R q . 

From these equations it can be seen that the upper level vortex is 

dependent on the increase in 0^ in the boundary layer. The problem 

becomes one of modeling the boundary layer air-sea interaction so that 

the 0^ increase corresponds to the upper level vortex. 
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B. CORGNAT I ' S MODEL 



1. Prediction Equation 

a. Atmospheric Prediction 

In the atmosphere the predicted variables were heat and 
moisture in the form of potential temperature (0^) and specific humid- 
ity ( c l a ) respectively. These two parameters, along with pressure (P), 

determined 9 which was the one factor in determining the maximum wind 
e 

(Equation 8). 

Partial derivatives were approximated using finite differ- 
encing techniques to solve the thermodynamic equation in the following 
form S0 [Sv 9 Sv R0~1 & 00 Q G 

St “i_w + iii + + c^p • (9) 

L - 1 P 

where cn is the vertical motion and Q is sensible heat transfer. The 

x s 

first term on the right hand side was neglected due to symmetry. The 

boundary layer was assumed well mixed with height, and therefore, 0^ 

could be calculated from the surface values of air temperature and 

pressure. Moisture changes were determined in a similar fashion using 

the expression Sq [Sv q Sv Rq 

a _ l e a . r a 

St [r§9 RoR 

where Q is the latent heat transfer, 
e 

Corgnati (1971) used bulk aerodynamic transfer formulas for 
determining the transfer of latent and sensible heat in the form 

% = f C p P s C D (T w - T a )V • (11) 

% - f L v P.Vv * ^a )V * (12) 

where is the specific heat of dry air, is the latent heat of 

vaporization of water, p g is the atmospheric density at the sea surface, 
T^ - T^ is the air-sea temperature difference, q^ and q^ are the 



Kl 



LAP 



( 10 ) 
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specific humidity of the air and the saturated specific humidity at the 
sea surface temperature. The surface wind (V) was approximated at 70$ 
of the wind at the top of the boundary layer (Equation 5). Furthermore, 
Robinson's (1966) heat transfer coefficient (2/3 C^) was used, which 
takes molecular processes near the surface into account. 

Vertical motion at the top of the boundary layer was com- 
puted from the continuity equation in polar coordinates. 




where the first term is zero due to symmetry. Vertical motion at the 
surface was assumed to be zero, and a boundary layer 100 mb thick was 
used, based on the assumption that effects within the boundary layer 
would not exceed this height. 

b. Oceanic Prediction 

Predicted variables in the oceanic portion of Corgnati's 
model were sea surface temperature (T^) and the saturated specific 
humidity at the sea surface temperature. A temperature gradient of 2C 
per 100 meters with a layer depth of 100 mete.rs described the ocean 
structure. The amount of heat removed, through latent and sensible 
heat transfer, was computed and the temperature of the ocean surface 
reduced accordingly. Corgnati found that, since the model was 
stationary, reductions in T^ led to reduced air-sea transfers and, 
therefore, reduced wind speed. In order to test only the atmospheric 
portion of the model the sea temperature was fixed. 

2. Initialization of Surface Temperatures and Humidities 

The relationship between 0 , tangential wind at the eye wall 

i 

(V ), and pressure at the eye wall (P.) made it possible, by specifying 

© . l 
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V , to solve for 9 and P. using Equations (7) and (8). The surface 

0. e t i 

pressure distribution was then determined by assuming a gradient wind 
balance. 

Initialization of the temperature and moisture fields, both 
in the boundary layer and at the ocean surface, were important parts 
of the model. The surface water temperature was set at 30C, then the 
available moisture at the ocean surface was calculated from the inte- 
grated form of Clapeyron’s equation 



Ln 



_s _ „ T-273 ? _ T 

6.105 25 T 5,2 L 273 



(14) 



where e g is the saturated vapor pressure. Air in contact with the 
water was assumed saturated, with the specific humidity calculated from 

.622 e 



P - .373e 



(15) 



The equivalent potential temperature was computed by applying the 
relationship between height gradient of the isobaric surfaces and the 
tangential wind field. 

By specifying the initial state from the equivalent potential 
temperature at the eye wall, the effect of expansion cooling and 
sensible heat source at the surface were taken into account. Riehl 
(1963) suggested that 9^ must rise above 350-352K to allow deep con- 
vection, which warms the column and lowers the central pressure to 
values which can sustain hurricane force winds. 

Specific humidity of the air (q^) was calculated assuming the 
initial air temperature (T ) equal to the water temperature, from 



q = T Ln , 

^a a L T ’ 

v a 

where T^ is the pseudo-equivalent air temperature. 



(16) 
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C. MODIFICATIONS TO CORGNATI’S NUMERICAL MODEL 



Corgnati (1971) used a 17 by 17 rectangular grid with 15 nautical 
mile spacing. The center of the grid was the hurricane center, and 
variables were computed along a diagonal and interpolated to the other 
grid points. Due to the symmetry of the model, the four quadrants had 
identical results, and as much information could be derived through the 
use of a one -dimensional model. Resolution was increased by using 100 
grid points spaced at 3 kilometer intervals along a radial. One dis- 
advantage of the increased resolution was that smaller time steps had 
to be used, although this was partially compensated by the reduction in 
the number of grid points. Inherent instability in a forward -in- time, 
centered-in-space finite differencing technique, made it necessary to 
change to a leapfrog technique with a forward time step every 24 time 
steps. 

The equivalent potential temperature gradient was the determining 
factor in predicting the maximum wind. When deriving Equation (8), 
Riehl assumed cyclostrophic flow, so that the relationship was only a 
good approximation in the area of maximum winds. Assuming gradient 
winds led to the relationship 



v •> < in 3 S9 

r + fv e 5 ST 



(17) 



which was used during the initialization to determine the equivalent 
potential temperature field. 

The interdependence of the equivalent potential temperature field 
and the wind field makes the selection of the eye wall and outer radii 
important factors in determining the wind field. The inner radius was 
allowed to adjust towards the radius of maximum equivalent potential 
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temperature, rather than remain fixed as in Corgnati's model, and the 
outer radius was determined from Equation (1). The maximum tangential 
wind component was then determined from Equation (17) using the gradient 
of equivalent potential temperature from the eye wall to the outer 
radius . 

Corgnati (1971) found that relative humidity values at the eye wall 
region exceeded saturated conditions after a short period of time. He 
theorized that cumulus clouds, in the form of "hot towers", with the 
moist ascending air being replaced by relatively dry subsiding air, 
would serve to remove this excessive moisture. To model the "hot towers" 
Corgnati removed all moisture in excess of a critical value of 95$ rela- 
tive humidity, and measured this as convective precipitation. A similar 
problem was encountered at the outer radius, where the computed value of 
0^ was found to be significantly higher than the 350K specified by Riehl. 
Since the maximum wind was dependent on the radial gradient of 0^, the 

increased 0 caused a decrease in the wind field. To correct this, a 
e 

critical relative humidity of 75$ was used at the outer radius. This 
simulated the relatively dry subsiding air outside a hurricane. Several 
profiles of critical relative humidity were experimented with and will 
be discussed later. 

Corgnati (1971) experimented with the release of latent heat, within 
the boundary layer, of the excess moisture removed through cumulus con- 
vection. Releasing all the latent heat within the boundary layer greatly 
increased the hurricane intensity. The bases of cumulus clouds within 
hurricanes have been observed to be lower towards the eye. To model 
this, latent heat release was made an inverse function of radial distance. 
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With these changes it was now possible, through latent heat re- 
lease, relative humidity control and ocean temperature adjustments, 
to model steady state hurricanes of different intensities. 



D. CARDONE'S BOUNDARY LAYER MODEL 



The use of bulk aerodynamic transfer formulas requires some know- 
ledge of heat and moisture exchange coefficients as well as winds 
within the surface layer. Recent advances in the study of atmospheric 
turbulence have led to a greater confidence in describing the distri- 
bution of wind and temperature in the atmospheric boundary layer. This 
progress has arisen largely through the application of the Monin-Obuhkov 
"similarity theory”, which has been supported by observational studies. 
Cardone (1969) has extended Blackadar's (1965) two-layer, neutral, 
baroclinic boundary layer model, over fixed terrain, to the boundary 
layer over the sea. ‘McConathy (1972) evaluated Cardone* s formulation 
of the marine boundary layer. The results of McConathy' s study were 
used to apply Cardone' s specification of the boundary layer, instead of 
the bulk aerodynamic transfer formulas in Corgnati's model, and to 
determine a radial wind profile throughout the boundary layer. 

The Monin-Obuhkov "similarity theory" for the constant stress, or 
surface layer, predicts that a universal relation should exist for the 
wind shear. 






kz u 



and the temperature gradient. 



(18) 



ae 

a 

<$z 



z *t 



where the friction velocity (U*) is defined as 

u * = $ 1/2 , 



(19) 



( 20 ) 
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and the scaling temperature 



®* ” kU, 



V. 



( 21 ) 



k is von Karman’s constant of .35, z is the height, H is the total heat 
flux, and ^ are functions of stability. Matching the integrated equa- 
tions for the Ekman layer to the integrated equation for the constant 
stress layer at top of the surface layer led to the following set of 
equations, in terms of and inflow angle (a), 

_2„ A ,h n i 1/2 



r - [2kB 0 sinZ<x i ( ^ )]J 

g 



(22) 



and 



k \/2 sin 



a) 



LnB Q Ro - t'Ch/L*) 



(23) 



where the integrated profile stability parameter is defined by 

l+<j) 1 1 i2 



1 - 6 - 3Ln{) + 2Ln(— — )+ 2 

T u r u 2 



.1 TT i+r 

tan <|) u -2 + Ln(-~i) 



(24) 



and the height of the surface layer (h) is defined by 

B V 

u - ° E 

h f » 



(25) 



the surface Rossby number (Lettau 1959) by 

V 

Ro = ft > 



(26) 



and the modified stability length (Panofsky 1963) by 






L = rr 1 L = 



U*T du/dz 
kG dQ/dz 



(27) 



In these expressions is the surface geostrophic wind, and B q is 
Blackadar's constant of 0.0003. The drag coefficient (C ) was defined 



in terms of the roughness parameter (Z ) and stability parameter as 

k ° - ! 



-[r 



z/Z Q - Uz/L 



- 7 ] 



(28) 
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The roughness parameter was formulated by Cardone (1969) as 



Z Q = 6,84 + 4.28 X 10' 3 U* - 4.43 x ltf 4 



( 29 ) 



where the second term on the right hand side represents Charnock' s (1955) 
formulation for aero dynamically rough flow under high wind speed condi- 
tions. 



.035in 



Z o = 



(30) 



Equations (22) and (23) were solved simultaneously for U^/V^ and Ot from 
the external parameters V , f, and air-sea virtual temperature differ- 



E. APPLICATION OF CARDONE' S BOUNDARY LAYER MODEL 

Cardone' s model was applied to the balanced hurricane to determine 
the fluxes of heat and moisture as well as momentum transfer. The 
determination of a mean radial wind component, using Cardone 's model, 
assured continuity in determining advection within the boundary layer. 

The tangential wind, which represented a gradient wind at the top 
of the boundary layer, was corrected for thermal wind and used as the 
surface gradient wind. 

To include the centripetal acceleration, the coriolis parameter was 
replaced by 

V 

f' = f + — . (31) 



By assuming the turbulent transfer coefficient for heat (K^) equals the 
turbulent transfer coefficient for momentum (K m ) , and knowing L'from 
the solution of Equations (22) and (23), the total heat transfer was 



determined using 



L * — L = . 
*h 



C p T 
P a 
kGH 



(32) 
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The total heat (H) was then separated into latent and sensible components 
by assuming their ratio was determined by 

q c /e - e \ 

_s _ _e (_§ « 



(33) 



'e “v ' ‘‘a 

The predictive Equations (9) and (10) required a radial wind for 
advecting the heat and moisture. The assumption that heat and moisture 
are homogeneous, within the depth of the boundary layer, made it possible 
to use an integrated value of radial wind. This was determined by 
solving the integrated form of Equation (18) 



” z =r r L ni- > 

o 



(34) 



at z = t i=l, 2, 3, 4, and averaging these values. The mean radial 
wind in the surface layer was then determined by multiplying by the sine 
of the inflow angle (a) . The inflow within the spiral layer was deter- 
mined through integration of the Ekman solution through the spiral layer 
to get co 

V V r dz= x /2 Vasina [cos(a+ 3 2 -Bh)-sin(a-l- 3 2 -Bh)],(35) 



where 



R . , ...f -,1/2 
B ( 2iych' 



(36) 



The two integrated values of from Equations (34) and (35) were then 
used in determining the vertical motion (Equation 13), and the advection 
of heat and moisture (Equations 9 and 10) in the predictive portion of 
the model. 
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III. DISCUSSION OF RESULTS 



The numerical process encompassed three stages. First was an ini- 
tialization phase in which the model was started with a given V^, T^, 

T , and R. , then the model went through an adjustment process in which 
a 1 

the moisture content in the boundary layer increased until either a 
balance was reached between latent heat flux and moisture advection, or 
the moisture reached the critical relative humidity values. The final 
stage was a "steady state" in which the time variations of heat and 
moisture in the boundary layer were small. 

The following discussion deals with the interpretation of results 
obtained when various parameters were changed, both in the initial and 
adjustment stages of the model. Due to the complex interactions of the 
model, a simplistic approach of establishing a base run and then varying 
one parameter at a time was used. This type of experimentation made the 
effects of the changes more obvious. 

A. CORGNATI'S MODIFIED MODEL 

The modifications to Corgnati's model, as discussed in the previous 
section, were incorporated to achieve more flexibility in describing 
hurricanes of varying intensity. Corgnati discussed in detail the adjust- 
ment process. He found that steady state was dependent on latent heat 
release and critical relative humidity used to simulate the cumulus con- 
vection. These parameters control, to a good part, the heat and 
moisture at the eye wall (R^), which in turn determines the equivalent 

potential temperature (9 ) and, therefore, the hurricane intensity 

i 

(Equation 8) . 
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The parameters varied were initial air temperature, water temperature, 
critical relative humidity, and latent heat release within the boundary 
layer. Figures 3 through 15 illustrate the relationships between the 
various parameters, and were obtained from the base run with water tem- 
perature set at 30C, air temperature initialized at 29C, critical rela- 
tive humidity set at 95$ at R^, with a linear profile to 75$ at R q , and 
no latent heat release within the boundary layer. 

Air temperature (Figure 3) exhibited a nearly isothermal (less than 
.3C) profile from R q to about 3R i? with an approximately logarithmic 
drop from 3R^ to a minimum temperature at R^. The temperature drop cor- 
responds to the' very strong horizontal pressure gradient (Figure 4) . 

The adiabatic cooling due to expansion would cause about a 5 C drop, 
but, as Riehl (1954) pointed out, the blowing spray in the high wind 
region would greatly increase the effective surface area and, therefore, 
the latent and sensible heat transfer. This enhanced transfer would off- 
set the adiabatic cooling and account for the slight temperature decrease 
towards the eye as observed by Deppermann (1939). The bulk transfer 
formulas (Equations 11 and 12) do not take into account the enhanced 
transfer due to the blowing spray. However, the increased wind speed 
and air-sea temperature difference were sufficient to increase the 
sensible and latent heat fluxes and offset much of the adiabatic cooling. 

Relative humidity (Figure 5) was at the critical value which indi- 
cates the steady state stage was one in which the flux and advection 
of moisture was kept in balance through the removal of moisture exceed- 
ing the critical relative humidity. Figure 15 indicates the moisture 
removed to keep the relative humidity at the critical value. 
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1 1 1 1 h 

60 120 180 240 300 

R km 

Figure 3. Radial profile of air temperature for the base run 

of Corgnati's modified model, T = 30C, T initialized at 29C. 

w a. 




Figure 4. Same as Figure 3 except pressure. 
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Figure 5. Same as Figure 3 except relative humidity. 




Figure 6. Same as Figure 3 except specific humidity. 
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Specific humidity (Figure 6) followed the relative humidity in the 
isothermal region. The sharp decrease in temperature reduced the 
amount of moisture, despite the increase in critical relative humidity, 
and the result was a peak in specific humidity outside the eye wall. 

Equivalent potential temperature (Figure 7) is dependent on tempera- 
ture, moisture, and pressure. In the nearly isothermal region and with 
small pressure gradient, 9^ is governed by the amount of moisture. From 
3R^ to the temperature decrease has a two-fold effect in that it also 
reduces specific humidity and 0^ becomes linear. The gradient of 9^ 

between R. and R determines V (Equation 17), and Equation (5) de- 

i 

scribes the tangential wind profile (Figure 8). Pressure and wind are 
linked through the gradient wind equation, and the sharp wind increase 
is reflected by the strong pressure gradient. 

In a hurricane the winds rotating about the eye experience both 
centripetal and coriolis acceleration (Figure 9) . In the outer portion 
of the hurricane coriolis dominates. However, as V increases and 
radius decreases, centrifugal force becomes dominant, supporting the 
strong pressure gradient in the eye wall regions, and causing a decrease 
in integrated V . The decrease in V q in the near surface region requires 
cross isobaric flow; this results in a frictional force opposing the 
pressure force. The importance of these various forces and their region 
of influence will become more apparent in the discussion of results 
obtained using Cardone's boundary layer model. Figure 10 shows the 
effect of cyclostrophic flow on the radial wind. The strong vertical 
motion (Figure 11) at the eye wall is due to the decrease in radial wind 
towards the eye wall and the decreasing circumference (Equation 13). 

In reality the vertical motion is due to deep convection maintained 
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Figure 7. Same as Figure 3 except equivalent potential temperature. 




Figure 8. Same as Figure 3 except tangential component of the wind. 
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Figure 9. Illustrating the forces acting in the boundary layer 
of an axisymmetric hurricane. 




Figure 10. Same as Figure 3 except radial component of the wind. 
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Figure 11, Same as Figure 3 except vertical motion. 
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through the buoyant energy gain represented by the 0^ increase from R q 
to R^. The convection is sufficient to cause the large pressure gradi- 
ent which is maintained by the cyclos trophic balance. 

Latent heat flux (Figure 12) is dependent on wind speed and the 
difference in specific humidity of the air and saturated air at the sea- 
surface temperature (Equation 12). The low moisture content of the air 

at R causes a local maximum in heat flux. The increase in moisture 
o 

dominates the region from R q to 3R^ causing a decrease in flux despite 
the increasing wind. Maximum heat flux occurs at R^ where the decrease 
in pressure and temperature acts to increase the specific humidity 
difference and the wind reaches a maximum velocity. 

Sensible heat flux (Figure 13), which is dependent on wind speed 
and air-sea temperature difference (Equation 11), reaches a maximum at 
R^ and is about an order of magnitude less than the latent heat flux. 

Precipitation was simulated by two mechanisms. One was the large 
scale moisture flux through the top of the boundary layer (Figure 14), 
which reflects the specific humidity and vertical motion. The second 
mechanism is the removal of moisture which exceeds the critical relative 
humidity (Figure 15) . Both curves are accumulated during the first 15 
hours and, therefore, are influenced by the adjustment phase. The peak 
in Figure 15 at R q is due to the movement of R q during adjustment. 

1 . Experiments with Temperature Variatio ns 

In a series of experiments the horizontally uniform water tem- 
perature was varied from 28 to 31C while air temperature was 
initialized at 28C. The radial profiles of the various parameters were 
similar to those obtained from the base run and illustrated in Figures 3 
to 15. There were, however, changes in magnitude and significant values 
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Figure 12. Same as Figure 3 except latent heat flux. 




Figure 13. Same as Figure 3 except sensible heat flux. 
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Figure 14. Same as Figure 3 except total moisture flux through the 
top of the boundary layer for first 15 hours. 




Figure 15. Same as Figure 3 except for total moisture removed by 
critical relative humidity mechanism for first 15 hours. 
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are presented in tabular form (Table 1) . Air temperature (T ) at R 

o 

adjusted to within one degree of the water temperature, and the air 

temperature (T ) at R. was about 1.5C cooler than T . One of the 
a i a o 

interesting relationships was the dependence of the model on the differ- 
ence in 9^ from R q to R^. Even when the water temperature was fairly 
cool, the increase in 0^ due to air-sea exchange was large enough to 

sustain hurricane force winds. Since 0 is representative of the ambi- 

o 

ent values in a tropical atmosphere, and values of 350K are necessary 

for deep convection, 0 was limited to a minimum value of 350K for 
o 

determining V (Equation 8). This limitation is reflected in the value 
°i 

of V for the 28C water temperature case. One of the assumptions of the 
i 

model is that the energy source for the hurricane is the ocean; there- 
fore, the hurricane should have more energy available and increase in 
intensity over warmer water. Intensity did increase from 49 m sec ^ for 
the 28C ocean to 62 m sec ^ for the 31C ocean. 

2. Experiments with Variations of Critical Relative Humidity 

Critical relative humidity was used to simulate decreases of 
moisture in the boundary layer due to convective mechanisms. To test 
the effect of the critical relative humidity, three profiles were speci- 
fied (Figure 16). Resulting relative humidity profiles (Figure 17) show 
that in the convex case relative humidity did not reach critical values 
except in the eye wall and outer regions. In the concave and linear 
cases relative humidity reached critical values at all radii between R q 
and R^. Latent heat flux is dependent on the difference in vapor 
pressure of saturated air at sea surface temperature and ambient air; 
because of this dependence the latent heat flux distribution was affected 
by the specification of critical relative humidity (Figure 18). 
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Table I 
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air temperature was initialized at 28C. 



T 
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T 

a. 

i 


T 

a 

o 


8 

e. 
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<°C) 


(°C) 


(°K) 


28 


26.7 


27.8 
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29 


27.0 


28.5 


363 


30 


27.6 


29.2 


367 


31 


28.4 


29.9 


371 



0 

e o 




P 

o 


(°K) 


(m/sec) 


(mbs) 


347 


49 


975 


350 


57 


965 


353 


59 


961 


356 


62 


958 
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Figure 16. Radial profiles of the three specified critical 
relative humidity curves. 




Figure 17. Radial profiles of relative humidity resulting from 
the three different specifications of critical relative humidity 
T w “ 30C, T a initialized at 29C. 
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Figure 18. Same as Figure 17 except latent heat flux. 




Figure 19. Same as Figure 17 except equivalent potential temperature. 
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Hurricane intensity, which is determined from 0 at R and R. , did not 

vary. The eye wall was made to adjust towards the maximum 9^; in the 

convex case maximum 0^ occurred outside the eye wall and R^ adjusted to 

36 km. In the concave and linear cases maximum 0 occurred at or inside 

e 

the eye wall and R^ adjusted to 24 km (Figure 19). 

It was concluded that the actual form of the critical relative 
humidity curve between fixed extrema was not crucial in predicting the 
intensity of the storm; however, it does affect the distribution of 
heat exchange. 

3. Experiments with Variations of Latent Heat Release 

Assuming dry adiabatic expansion in the inflow would result in 
air temperature in the eye wall region as much as 5C colder than ob- 
served temperatures. Riehl (1954) hypothesized that latent heat release 
or sensible heat flux, enhanced by blowing spray, must offset the expan- 
sion cooling. To simulate the effects of this enhanced heating, the 
latent heat of condensation of the moisture released to maintain relative 
humidity below the critical value (Figure 15), was added to the boundary 
layer. Since the enhanced heating would be a function of the amount of 
blowing spray, and the amount of blowing spray would be dependent on 
wind speed, and wind speed is a function of radius, the percentage of 
heat released was made a linear function of radius with the maximum value 
at R^, and no latent heat release at R q . Corgnati experimented with re- 
leasing 100$ of this latent heat at all radii and found a marked increase 
in intensity. Table II illustrates the effects of latent heat release in 

the modified model. The increase in 0 , while 0 remained unaffected, 

e. e 

l o 

resulted in an increased intensity. 



38 



Table II 



Variations of significant values with latent 
heat released in the boundary layer. 



latent heat 
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(°C) 
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58 


963. 


25 


28.3 
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951 


50 


28.5 


29.5 
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75 


28.7 


29.5 


375 


354 


72 


941 


100 


28.6 


29.6 


376 


354 


74 
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Although latent heat was used, the heat gain could also be pro- 
duced by an enhanced sensible heat transfer due to blowing spray. In 
all probability both latent heat release and enhanced sensible heat 
transfer play some role in offsetting expansion cooling in the boundary 
layer. The temperature profile of the model, without latent heat re- 
lease, was in agreement with observed temperature profiles. The model 
was therefore run with no latent heat release in the boundary layer. 

B. CARDONE'S BOUNDARY LAYER MODEL 

Cardone's two layer baroclinic boundary layer model assumes hori- 
zontal homogeneity, stationarity, equal to a constant in the Ekman 
layer and 

K m = kU*z/ 4> u ( Z /L' ) (37) 

in the surface layer. The model also involves various empirically 
determined constants based on data from mid-latitudes. Despite these 
limitations it was felt that a boundary layer based on similarity theory 
would provide an internally consistent model of an hurricane boundary 
layer. The inputs to Cardone's model were surface gradient wind, 
virtual temperature difference between the air in the boundary layer and 
saturated air at the sea surface temperature, and latitude. 

1 . Application of Cardone's Boundary Layer Mode l 

Several problems arose with the direct application of Cardone's 
model in place of the bulk transfer equations and the integrated inflow 
from the tangential wind equation. The model hurricane decayed rapidly 

during the adjustment phase, and heat and moisture fluxes were well 

-2 -1 

below the 3000 to 5000 cal cm day values determined by Whitaker 
(1967) to occur in the hurricane force wind region of hurricane Betsy. 
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Investigation showed that radial wind (Figure 20), determined from Equa- 
tions (34) and (35), increased inward from 3R^ to EL. Because V^_ must 
go to zero in the eye, vertical motion at the eye wall was on the order 
of 1500 mb hr ^ . In this region the vertical transport of moisture ex- 
ceeded the sum of the horizontal transport and evaporation, which de- 
creased 0 and led to a decrease in intensity. Haltiner and Martin 
i 

(1957) have shown that the balance of forces in a vortex may be depend- 
ent on centrifugal force. This was included in Cardone's model by 
replacing f by f f as formulated in Equation (31). The time scaling 
factor then became 

(f + V^' 1 )" 1 , (38) 

which had the effect of reducing the height of frictional influence as 
V q R ^ increased. In particular, the height of the constant stress layer 
was reduced from about 50 meters at R q to about 20 meters at R^ (Figure 
21). To maintain consistency within the boundary layer, the centripetal 
acceleration also had to be included in term B (Equation 36) of the Ekman 
spiral. The formulation for within the spiral layer became 




where V q (z) is defined by 




To determine in the spiral layer an iterative technique was employed 
in determining V 0 (z) using heights (z) at 10 meter intervals. Then V r (z) 
was determined at each height (Figure 22) using the V Q (z) for that height. 
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Figure 20. Radial profile of radial component of the wind resulting 
from direct application of Cardone's boundary layer model. 




Figure 21. Radial profile of height of the surface layer, illustrat- 
ing decreasing height of frictional influence due to cyclostrophic 
flow. 
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Figure 22. Vertical variation of radial component of the wind in 
the spiral layer, derived using iterative technique. 
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The integrated radial wind in the spiral layer was then determined by 
averaging the values of V^(z) over the spiral layer. Using an iterative 
technique proved to be rather time consuming. Because V q (z) increased 
rapidly with height, the assumption that V q (z) was equal to at the top 
of the boundary layer resulted in integrated values of V^(z) that closely 
approximated the averaged values of V^(z) (Table 3). The assumption 
that V q was constant with height in the spiral layer was used in deter- 
mining the integrated value of (Figure 23). Including centrifugal 
forces was now consistent within the entire model and the rapid decay 
described earlier, as well as the extreme vertical motion, was no longer 
evident. 

In a study of water temperatures in the wake of hurricane Betsy, 

Whitaker (1967) estimated the heat flux in the hurricane force wind 

-2 -1 

region to be between .035 and .058 cal cm sec . The application of 

Cardone's model produced heat fluxes with a maximum value around .02 cal 
-2 -1 

cm sec (Figure 24). One of the assumptions in Cardone's model was 

that the ratio of the turbulent transfer coefficient for heat (K^) to 

turbulent transfer coefficient for momentum (K^) was unity. Businger et 

al (1971) found ”The ratio of the eddy dif fusivities (a^) for heat and 

moisture is greater than unity. The value at neutral is 1.35; and shows 

a marked increase with instability.” Since heat transfer, as determined 

in the model, is directly dependent on a, (Equation 32), increasing a, 

n h 

increased the heat flux. In a hurricane the large quantity of blowing 
spray should serve to enhance the heat transfer. Furthermore, the 
velocity imparted to this blowing spray should serve to increase the 
momentum transfer. To simulate the increasing heat flux and momentum 
transfer, the roughness length (Z Q ) was increased by increasing the 
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Table III 



Comparison of integrated V values obtained 
using iterative technique, and assuming V (z) 
constant. 



Radius (km) 

V (m/sec) 
iterative 
technique 

V (m/sec) 

V (z) const 
© 



30 


60 


90 


120 


150 


180 


210 


3.14 


3.22 


3.42 


3.49 


3.40 


3.38 


3.20 



3.22 3.27 3.47 3.33 3.43 3.40 3.21 




Figure 23. Radial profile of integrated radial component of the 
wind assuming cyclostrophic flow and V constant within the spiral 
layer. 0 
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Figure 24. Radial profile of combined sensible and latent heat flux 
obtained using Cardone's model, C« = .00428, a, = 1.0, T = 30C, 

and T initialized at 29C. 
a 




Figure 25. Same as Figure 24 except a h = 2.0, = .05. 



46 



2 

the coefficient (C.) of the U. term in Cardone's formulation of Z (Equa- 
Z * o 

tion 29). It was found that realistic values of heat flux were obtained 
when C ^ was increased from .00428 to .05, and a^ taken as 2.0 (Figure 25). 

Cardone's model, with cyclostrophic flow, C ^ increased to .05, 
and a^ increased to 2, was now used to determine the sensible and latent 
heat fluxes, replacing Equations (11) and (12), and the average radial 
wind component, replacing Equation (6). 

A base run was made in which water temperature was set at 30C, 
air temperature initialized at 29C, critical relative humidity set at 
95$ at with a linear decrease to 75$ at R q , and no latent heat re- 
leased within the boundary layer. This base run was compared to the 
base run illustrated by Figures 3 through 15 and discussed earlier. Car- 
done's model produced greater latent and sensible heat flux (Figures 12 
and 25) which led to slightly warmer temperatures at R^ and increased 

0 . As a result of the increased 0 , V was about 5 m sec ^ greater, 

e. e. o. & 

i ii 

Table 4 compares some of the significant values of the base runs. The 
use of Cardone's boundary layer model provided some additional para- 
meters that were indicative of the structure of the boundary layer, as 
illustrated in Figures 26 through 29. The modified stability length 
(L 1 ), illustrated in Figure 26, can be interpreted to be proportional to 
the ratio of mechanical and buoyant production of turbulent kinetic 
energy. The large increase of L', with increasing V 0 , reflects an in- 
crease in the mechanical versus the buoyant production of kinetic energy. 
Surface stress values (Figure 27), which increases rather dramatically to 
a maximum value at R^, are in close agreement with stress values deduced 
by Miller (1962) using data from hurricane Helene (Table 5). The increase 
in stress towards the center is due to the increase in wind and the 
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Table IV 



Comparison of significant values obtained from 
two base runs. 

Corgnat i ' s Cardone ' s 



parameter 


modified model 


model 


V (m/sec) 

i 


58 


62 


T (°C) 

a i 


27.9 


28.0 


T (°C) 

a 

o 


29.5 


29.2 


e - e (°k) 

e. e 

1 o 


13 


15 


P (mbs) 


963 


953 



48 



— 1 500 — 



— 1000 
L' 

mefers 



- 500- 




R: 60 120 , 180 240 R 

K km 

Figure 26. Same as Figure 25 except modified stability length. 




Figure 27. Same as Figure 25 except surface stress. 
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Table V 



Comparison of surface stress determined 
by Miller, and values obtained from the 
hurricane model. 



surface stress 


Radius 


0-20nm 


20-40nm 


40-60nm 


dynes /cm^ 


Miller's 

values 


60.0 


51.0 


27.6 




Radius 


20nm 


30nm 


50nm 




computed 


60 


35 


20 



50 



rougher surface. The increase of stress with decreasing radial distance 
is reflected in the increase in inflow angle within the constant stress 
layer (Figure 28); also in the radial velocity at the top of the sur- 
face layer (Figure 29). These values are comparable to inflow angles 
and radial velocities observed in actual hurricanes. One of the assump- 
tions made in Corgnati's model was that the cosine of the inflow angle 
could be closely approximated by unity. Despite the increase in the 
inflow angle, this appears to be a good approximation. 

2. Temperature Experiments with Cardone's Model 

To further test the use of Cardone's boundary layer model, ex- 
periments were conducted in which water temperature was varied from 28 
to 31C. These experiments were similar to the ones conducted with 
Corgnati’s modified model, and discussed in section IIIA1. The heat 
flux distribution obtained, using Cardone's model, appeared to be more 
sensitive to temperature variations than the heat flux distribution 
obtained using the bulk aerodynamic transfer formulas (Figure 30). Due 
to the differing heat flux, the variation in the difference between 
extrema values of 0^ was greater, using Cardone's model, than the 
variations obtained using bulk aerodynamic transfer formulas. The 
result of these greater variations was that, using Cardone's model, a 

sea surface temperature increase from 28 to 31C produced a V increase 

-1 0i 
of 25 m sec (Table 6); while using bulk aerodynamic transfer formulas, 

the same temperature increase produced only a 13 m sec ^ increase in 

V & • It was felt that the greater sensitivity to sea surface tempera- 
i 

ture, exhibited by Cardone's boundary layer model, more closely simulated 
the dependence of a hurricane on the oceanic heat source. 
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Figure 28. Same as Figure 25 except inflow angle in surface layer. 




Figure 29. Same as Figure 25 except radial component of the wind 
at top of surface layer. 
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Figure 30. Radial profiles of combined latent and sensible heat 
flux, comparing values obtained from Cardone's model with T = 31C 
and T = 29C, and Corgnati's modified model with T = 31C and 

T = W 29C, T initialized at 28C. w 

w a 



Table 



Changes in significant 
ing water temperature, 
initialized at 28C. 
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IV. CONCLUSIONS 



The results of this study have indicated that a marine boundary 
layer based on similarity theory can be applied to an hurricane and 
improve the effects of air-sea interaction. 

First, results indicate that the time scaling factor in cyclo- 
s trophic flow is dependent on (f + V q R \ rather than f as in 
Cardone's model. The effect of the centripetal acceleration was to 
reduce the height of frictional influence as wind speed increased and 
radial distance decreased. 

Second, observed temperature fields indicate that enhanced heat 
flux is necessary to offset the adiabatic expansion cooling within 
the boundary layer of an hurricane. This enhanced heat flux is likely 
to be due to the large quantities of blowing spray observed in an 
hurricane. Consequently, the ratio of the coefficients of heat dif- 
fusivity to momentum diffusivity was increased from unity to 2. 
Increasing the roughness length, as a function of friction velocity 
squared, was used to simulate the increased momentum transfer and 
resulted in surface stress values that appear to be realistic. The 
increase in roughness length was further justified by resultant heat 
flux values that were in agreement with determined values. 

Third, the removal of moisture from the boundary layer, through 
convective mechanisms, appears to have a strong influence on air-sea 
interaction. Vertical advection by the mean motion is not sufficient 
to export the large horizontal flux of moisture within the boundary 
layer and the evaporation. The removal of moisture in excess of a 
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critical relative humidity value appears to adequately simulate f, hot 
tower" convection and surrounding subsidence. 

Finally, in the earlier model of Corgnati (1971) it had been neces- 
sary to specify a drag coefficient and extrapolate a surface wind to 
determine air-sea interaction. However, the external parameters 
required to use similarity theory were surface gradient wind, and the 
virtual temperature difference between saturated air at the sea surface 
temperature, and the air in the boundary layer. The vertical structure 
of the wind and an implied drag coefficient were inherent in the 
iterative solution based on similarity theory. 

The more complete boundary layer specification based on similarity 
theory appears to be justified in view of the good agreement obtained 
between predicted and observed or estimated values of radial inflow, 
heat flux, and surface stress. Furthermore, the boundary layer model 
calculates wind profiles within the surface layer, as well as surface 
stress; and effects of both stability and changes of surface roughness 
on air-sea interaction have been included in the hurricane model. 

One of the assumptions made in this study was that the stability 
length was defined by density stratification based on virtual tempera- 
ture. This assumption was crucial in the determination of latent heat 
flux. 

Further studies, in the area of air- sea interaction under high wind 
speed conditions, are necessary to provide guidance in the selection of 
empirical constants in the similarity solution. Specifically, the 
effects of blowing spray, in enhancing heat flux and momentum transfer, 
may be important in maintaining tropical cyclones. 



55 



BIBLIOGRAPHY 



1. Blackadar, A. K. , 1965: A simplified two-layer model of the baro- 

clinic neutral atmosphere boundary layer. Air Force Cambridge 
Research Laboratories Report 65-531, pp. 49-65. 

2. Businger, J. P., J. C. Wyngaard, U. Izumi, and E. E. Bradley, 

1971: Flux-profile relationships in the atmospheric surface 

layer. J. Atmos. Sci ., 28, 181-189. 

3. Cardone, V. J., 1969: Specification of the wind distribution in 

the marine boundary layer for wave forecasting. New York 
University, School of Engineering and Science, Scientific Report 
GSL-TR69-1, University Heights, New York. 

4. Chamock, H., 1955: Wind stress on a water surface. Quart . 

J, Roy. Meteor. Soc ., 81, 639. 

5. Corgnati, L. B., 1971: Experiments with a simple hurricane- 

interacting ocean model, M.S. Thesis, Naval Postgraduate School, 
Monterey. 

6. Deppermann, C. E., 1939: Some characteristics of Philippine 

typhoons. Manila, Bureau of Printing, 143. 

7. Haltiner, G. J. , F. L. Martin, 1957: Dynamical and Physical 

Meteorology . New York, McGraw-Hill Book Co., 181-183, 234-241. 

8. Lettau, H. H., 1959: Wind profile, surface stress and geostrophic 

drag coefficients in the atmospheric surface layer. Advances in 
Geophysics , 6, 241-257. 

9. McConathy, D. R. , 1972: Specification of marine surface boundary 

layer parameters, M.S. Thesis, Naval Postgraduate School, Monterey. 

10. Miller, B. I., 1962: On the momentum and energy balance of 

hurricane Helene (1958). National Hurricane Research Project 
Report No. 53, 10-13. 

11. Panofsky, H. A., 1963: Determination of stress from wind and 

temperature measurements. Quart. J. Roy. Meteor. Soc ., 88, 85-94. 

12. Riehl, H. , 1954: Tropical Meteorology . New York, McGraw-Hill 

Book Co., 281-356. 

13. Riehl, H. , and J. S. Malkus, 1961: Some aspects of hurricane 

"Daisy", 1958. Tellus 13, 181, 213. 

14. Riehl, H. , 1963: Some relations between wind and thermal structure 

of steady-state hurricanes. J. Atmos. Sci ., 20, 276-287. 



56 



15. Robinson, G. D. , 1966: Another look at some problems of the air- 

sea interface. Quart. J. Roy. Meteor. Soc ., 92, 451-465. 

16. Shea, D. J. , 1972: The structure and dynamics of the hurricanes 

inner core region. Department of Atmospheric Science, Colorado 
State University, Atmospheric Science Paper No. 182, Fort Collins, 
134. 

17. Whitaker, W. D., 1967: Quantitative determination of heat trans- 

fer from sea to air during passage of hurricane Betsy, M.S. 

Thesis, Graduate College of the Texas A&M University, 65, 28. 



57 



INITIAL DISTRIBUTION LIST 



No. 



1. Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 

2. Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 

3. Professor R. L. Elsberry, Code 51Es 
Department of Meteorology 

Naval Postgraduate School 
Monterey, California 93940 

4. Professor D. F. Leipper 
Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 

5. Lieutenant Commander N. A. S. Pearson 
Helicopter Support Squadron Seven 
Naval Air Station Imperial Beach 
Imperial Beach, California 92032 

6. Naval Weather Service Command 

Naval Weather Service Command Headquarters 
Washington Navy Yard 
Washington, D. C. 20390 

7. Lieutenant Commander R. W. Lyons 
Fleet Weather Central 
Norfolk, Virginia 23511 

8. Professor K. L. Davidson, Code 5lDs 
Department of Meteorology 

Naval Postgraduate School 
Monterey, California 93940 

9. Professor R. L. Haney, Code 5lHy 
Department of Meteorology 
Naval Postgraduate School 
Monterey, California 93940 

10. Professor R. J. Renard, Code 51Rd 
Department of Meteorology 
Naval Postgraduate School 
Monterey, California 93940 



Copies 

2 

2 

8 

1 

2 

1 

1 

1 

1 

1 



58 



Security Classification 



DOCUMENT CONTROL DATA - R & D 

(Security classification of title, body of abstract and indexing annotation must be entered when the overall report Is classified) 



t originating ACTIVITY (Corporate author) 


2a, REPORT SECURITY CLASSIFICATION 


Naval Postgraduate School 


Unclassified 


Monterey, California 93940 


26. GROUP 


3 REPOR T TITLE 



Further Experiments with an Hurricane-Ocean Interaction Model 



4. DESCRIPTIVE NOTES (Type of report and. Indue i ve dates) 

Master's Thesis; September 1972 

5- AuThORISI ( First name, middle initial, lest name) 

Nils A. S. Pearson 



6 REPORT DATE 

September 1972 


7 a. TOTAL NO. OF PAGES 
60 


76. NO. OF REFS 

17 


14. CONTRACT OR GRANT NO. 


9*. ORIGINATOR'S REPORT NUMB E R(S) 


6. PROJECT NO- 






C. ' 


96. OTHER REPORT NOI3I (Any other numbara that may be aa a I gned 
thia report) 


d. 




- 



10 DISTRIBUTION STATEMENT 



Approved for public release; distribution unlimited. 



1 1. SUPPLEMENTARY 



13 ABSTR 



NOTES 



12. SPONSORING MILITAR' 



ACTIVITY 



Naval Postgraduate' School 
Monterey, California 93940 



Modifications were made to Corgnati's (1971) steady state, symmetrical 
hurricane model, which was based on a model proposed by Riehl (1963), and 
in which bulk aerodynamic transfer formulas were used to predict the air- 
sea interaction below a specified vortex. The modifications resulted in a 
more flexible model in which hurricane intensity was controlled by sensible 
and latent heat fluxes, which were in turn dependent on the ocean surface 
temperature. A boundary layer model based on similarity theory was then 
used, to specify air-sea interaction, by replacing the bulk transfer formulas 
with Cardone' s (1969) extension of Blackadar's (1963) two-layer neutral, 
baroclinic boundary layer model. To achieve realistic results the time scale 
was modified to include centripetal acceleration, which resulted in a decreas- 
ing height of frictional influence with decreasing radius. Varying Cardone' s 
formulation of roughness length and the ratio of eddy exchange coefficients 
for heat and moisture resulted in stress and heat flux values in agreement 
with values determined from hurricane observations. Including the complete 
boundary layer model resulted in greater sensitivity to ocean temperature 
variations than when bulk transfer formulas were used to determine the air- 
sea interaction. 



DD,T;..1473 lPAGE '» 

S/N 01 01 -807-681 1 59 



Security Classification 



Security Classification 



t 4 


KEY WORDS 


LIN 


K A 


LIN 


K B 


L IN 


K C 






ROLE 


W T 


ROLE 


W T 


ROLE 


W T 


Air-sea interaction 
Hurricane model 
Boundary layer 
Latent heat release 
Sensible heat flux 















DD , F ho R v m „1473 <back) 

S/N 01 o t • 807- 6 8 2 1 60 Security Classification 




























I 



i 3 SEP 73 



2 1 6 u 0 



Thesis 138060 

P316 Pearson 

c.l Further experiments 

with an hurricane-ocean 
i nteraction model . 

• •• SEP 73 2 1 6 u n 





Thesis .UC<0b0 

P316 Pearson 

c.l Further experiments 

with an hurri cane-ocean 
i nteraction model . 



